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TRANSONIC FLUTTER INVESTIGATION OF TWO 64° DELTA 
WINGS WITH SIMULATED STREAMWISE RIB AND 
ORTHOGONAL SPAR CONSTRUCTION 
By George W. Jones, Jr., and Lou S. Young, Jr. 

SUMMARY 


An experimental investigation has been made in the Langley transonic 
blowdown tunnel of the transonic flutter characteristics of two 64° swept- 
back delta wings. Each wing had simulated streamwise ribs and orthogonal 
spanwise spars, but one wing had a different stiffness and mass than the 
other. Flutter was obtained on the more flexible wing at several Mach 
numbers from 0.79 to 1.28 and on the stiffer wing at Mach numbers from 
0.84 to 0.97. 

At a given Mach number, the value of the mass ratio at flutter dif- 
fered for the two wings by a factor up to 2, but the data were correlated 
by use of a parameter consisting of the flutter-speed coefficient divided 
by the square root of the mass ratio. As the Mach number was increased, 
the dynamic pressure required for flutter and the flutter frequency 
increased by a factor of about 2 at a Mach number of approximately 1.05) 
this increase is interpreted as a change in flutter mode. Reference 
flutter speeds were calculated by use of streamwise two-dimensional 
incompressible aerodynamic coefficients in a coupled modal analysis. 

These calculated flutter speeds were too high for the low-frequency 
flutter mode by 11 to 28 percent and too low for the high-frequency 
flutter mode by 11 to 35 percent. 


INTRODUCTION 


Although delta-wing plan forms are in present-day and projected use 
for high-speed aircraft, little is known about the transonic flutter char- 
acteristics of this type of wing configuration. An exploratory investi- 
gation has accordingly been made in an attempt to define some of the tran- 
sonic flutter problems of delta wings. The investigation, which was made 
in the Langley transonic blowdown tunnel, consisted of transonic flutter 
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tests on a 64° delta-wing plan form (aspect ratio about 2) which simu- 
lated in a crude manner one of several general types of delta-wing 
construction in present-day use. 

The model construction simulated spars normal to the fuselage plane 
of symmetry and streamwise ribs. Flutter points were obtained in the Mach 
number range from 0.79 to 1.28 with the model cantilever mounted at zero 
angle of attack without body freedoms. The effect of variations in the 
mass ratio was determined at subsonic speeds by use of a second, stiffer 
model which also had simulated streamwise ribs with orthogonal spars. 

The results of the investigation are presented herein together with 
a comparison of the experimental flutter speeds with those calculated by 
a simplified method. Also included is a comparison of the measured vibra- 
tion modes with those calculated from measured structural influence 
coefficients . 


SYMBOIS 


a streamwise distance from strip reference axis to strip 

center-of -gravity location, positive if center of gravity 
is behind reference axis, ft 

b streamwise strip semichord passing through influence coef- 

ficient stations on strip, ft 

t) r reference wing streamwise semi chord, mean geometric exposed 

semi chord, ft 

measured coupled natural frequencies (i = 1, 2, 5, or 4), cps 

g-^ structural damping coefficient in bending 

h. contribution of the ith linearized mode to nondimens ional 

vertical displacement of a wing strip reference axis 
(i = 1, 2, or 3) 

I a S mass moment of inertia of streamwise wing strip of width 8 

about the strip reference axis , slug-ft2 


M Mach number 

m 5 mass of a streamwise wing strip of width 8, slugs 

mass of wing sections (s = 1, 2, 3> • • •> 20), slugs, (fig. 4) 
s 

q dynamic pressure, ^ pV 2 , lb/sq ft 
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Sa5 static mass moment of streamvd.se wing strip of width 5 about 

strip reference axis, positive if center of gravity is 
behind reference axis, slug-ft 

t/c ratio of wing thickness to streamwise chord 

V free-stream velocity, ft/sec 

Ve/VR flutter-speed ratio, ratio of experimental flutter speed to 
calculated, or reference, flutter speed 

x Q g streamwise distance from leading edge to center-of-gravity 

position, fraction of streamwise chord 

y vertical deflection of wing 


y s 


a i 


5 

n 


vertical deflection of wing at a point where an influence 
coefficient was measured (s = 1 , 2 , . . . , 20 ) 

contribution of the ith linearized mode to torsional deflection 
of a wing strip about the strip reference axis (i = 1 , 2 , or 5 ) 

width of wing strip used for reference flutter-speed calcula-r 
tions, ft 

nondimens ional distance along exposed wing span, 

Spanwise distance measured from wing root 
Length of exposed span 


M- 

A 

P 

oo 

o>l 


mass-ratio parameter, 


Exposed panel mass 

(Exposed panel span)^Ttpb r ^] 


sweepback angle of leading edge, deg 

air density, slugs/cu ft 

angular frequency of flutter, radians/sec 

angular coupled natural frequencies, 2 jtf^, radians/sec 

angular coupled predominately torsion frequency, 005 , 
radians/sec 
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Subscripts: 

e experimental values 

R calculated values 


MODEIS 


Wing Geometry 

A sketch of both delta wings showing plan forms, basic dimensions, 
and construction is given in figure 1. Each wing had a leading-edge 
sweepback angle of approximately 64° and tips clipped along streamwise 
lines. The wing sections had a rounded leading edge over the first 
4 percent chord, straight parallel top and bottom surfaces to 85 percent 
chord, and a straight taper on top and bottom from 85 percent chord to 
a sharp trailing edge. Along the span each wing panel had a nearly con- 
stant ratio of thickness to chord except that near the tip the thickness 
ratio increased somewhat. (See fig. 2.) 


Wing Construction 

Each of the wings was constructed from a blank of 2024 aluminum 
alloy which was shaped into two panels, as described in the previous 
section, with an integral mounting block, as shown in figure 1. One 
panel of each wing was modified so as to roughly simulate spars normal 
to the fuselage plane of symmetry and streamwise ribs. The formation 
of the simulated ribs and spars was accomplished by cutting a pattern 
of circular holes, some of which were connected by streamwise saw cuts, 
through the solid 2024 wing panel. On both modified wing panels the 
holes and cuts were filled with lightweight, low— stiffness foam plastic 
and wrapped with a sheet of 0.005 -inch-thick fiber glass. Figure 1 shows 
that wing 2 has a pattern of streamwise saw cuts which is different from 
that of wing 1. The difference in patterns of the saw cuts changed the 
number and location of the simulated spars and, along with a slightly 
greater thickness for wing 2 , resulted in wing 2 being stiff er than wing 1 . 


Wing Physical Parameters 

Measurements were made on wings 1 and 2 of the first four coupled 
natural frequencies and node lines, the exposed panel mass, and the 
structural damping coefficient in bending. Values of the structural 
damping coefficient were determined from the decrement of free-bending 
vibrations in still air. The measured frequencies and node lines which 
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are presented in figure 3 were obtained by use of an electromagnetic 
shaker mounted close to the root (see fig. 3)j sprinkled salt crystals 
depicted the node lines at the natural frequencies. On the basis of 
repeatability of the data the measured frequencies have an accuracy 
which varies from about 1 percent for the lowest frequency to about 
2 percent for the fourth frequency. The calculated node lines and 
frequencies of wing 1 which are shown in figure 3 will be discussed 
sub s e quent ly . 

Structural influence coefficients were measured at 15 points on 
wing 1. The values of the influence coefficients obtained are given 
in table l(a). The points at which the influence coefficients were 
measured are shown in figure k which will be discussed in more detail 
later . 

Each influence coefficient was obtained in the following manner: 

The wing was firmly clamped in a horizontal plane to a massive mount. 

A traveling overhead support held a depth micrometer which was used 
to measure the deflections of the wing. The micrometer could be read 
directly to the nearest ten-thousandth of an inch and interpolated to 
the nearest fifty-thousandth of an inch. In order to ascertain when 
the pointed end of the micrometer (radius about 0.015 inch) touched 
the wing surface, a direct current electrical circuit containing a 
neon test lamp was rigged between the wing surface (which was coated 
with a conducting silver paint) and the micrometer point. When the 
point touched the wing, the neon test lamp lit and a reading was taken. 

A spark jump as the point neared the wing was virtually eliminated by 
a high resistance in the circuit. A micrometer reading was made of the 
wing position with no load on the wing, then a weight was hung at the 
desired station and another reading was taken. The loading-deflection 
reading procedure was repeated several times for each influence coef- 
ficient. The deflections thus obtained were averaged and then adjusted 
to give deflection per unit load. 

For use in frequency calculations the matrix of table l(a) was made 
symmetrical by the following procedure: Each influence coefficient of 

table l(a) was weighted (multiplied) by the number of deflections averaged 
to obtain it. Then the sum of each weighted pair of supposedly reciprocal 
influence coefficients was divided by the total of the weighting factors 
to give the final influence coefficient value. Influence coefficients at 
five additional stations (fig. 4) near the root were interpolated from the 
adjusted measured values and assumed zero deflection at the root. The final 
matrix of adjusted and interpolated influence coefficients used for fre- 
quency calculations is given in table l(b). In table l(a), 'jQ percent 
of the off-diagonal elements were within 3 percent of their corresponding 
adjusted values in table l(b), and 92 percent were within 6 percent. 
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After testing, the exposed modified panel of wing 1 was cut into six 
streamwise strips; and the mass, the mass moment of inertia about the 
assumed strip reference axis, and the center-of -gravity location of each 
strip were determined. The methods used to measure these parameters are 
discussed in reference 1. Each of the wing strips was then cut into 
sections associated with the influence coefficients and the masses and 
center-of -gravity locations of the sections were determined. 

The division of wing 1 into strips and sections, the points at which 
influence coefficients were measured and inferred, the assumed strip refer- 
ence axes, the section center-of -gravity locations, and a list of section 
masses are given in figure 4. Table II gives physical parameters of the 
wings as follows: Table Il(a) lists some basic physical properties of 

wings 1 and 2; table Il(b) lists the measured mass properties of wing lj 
and table Il(c) lists some computed deflection properties of wing 1 which 
will be discussed subsequently. 


APPARATUS AND TESTS 


The instrumentation, tunnel characteristics, and testing technique 
are described in detail in reference 1; only a brief description of these 
items is given in the following paragraphs. 

The flutter tests were made in the Langley transonic blowdown tunnel 
which has a slotted octagonal test section measuring 26 inches between 
flats. During the operation of the tunnel, a selected Mach number from 
subsonic Mach numbers up to supersonic Mach numbers of about 1.4, which 
is set by an orifice plate downstream of the test section, can be held 
approximately constant (after the orifice is choked) while test-section 
pressure, and thus density, is varied. The density range is approximately 
0.001 to 0.012 slug per cubic foot. 

The delta-wing models were cantilever mounted at 0° angle of attack 
in a cylindrical sting fuselage mount which covers the mounting block 
shown in figure 1. The sting fuselage mount extends without change of 
diameter into the subsonic-flow region of the tunnel and thus prevents 
the formation of bow shock waves which might reflect from the walls onto 
the model. The fundamental bending frequency of the sting fuselage mount 
with model attached is approximately 15 cycles per second. 

Basically, the instrumentation was as follows: Wire strain gages, 

located as shown in figure 5, were used to indicate model deflection about 
two different axes. A recording oscillograph was used to obtain contin- 
uous records of the strain-gage signals, tunnel stagnation temperature and 
pressure, and test-section static pressure. The records of the strain- 
gage signals were used to determine the start of flutter and the frequency 
of wing oscillations. 
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The Mach number range over which flutter was obtained on the wings 
was from 0.79 to 1.28, but an attempt was made to flutter the wings at 
Mach numbers up to 1.35 and dynamic pressures up to 4,500 lb/sq fb. 


ANALYSIS 


Calculation of Frequencies and Mode Shapes 

The symmetrical matrix, table 1(b), of adjusted and interpolated 
influence coefficients measured on wing 1 was put into the matrix equation: 

( y s j = <a 2 M[“s] {^s} (1) 

where I is the matrix of influence coefficients. Equation (l) was 
solved for the frequencies of the first three coupled natural vibration 
modes and for the nondimensional vertical deflections of each of the 
20 influence-coefficient stations of wing 1 in these three modes. As 
shown in figure 4, the 20 influence-coefficient stations of wing 1 were 
generally not located exactly at the center-of -gravity positions of the 
influence-coefficient sections. As a check, a matrix of influence coef- 
ficients at the center-of -gravity positions of the wing sections was 
inferred graphically from the values of table l(b) and is presented in 
table III. The matrix equation (eq. (l)) using this matrix was solved 
for the frequencies and wing section deflections of the first four natural 
vibration modes. A comparison of the measured frequencies with the two 
sets of calculated frequencies follows: 


Mode 

Measured 

frequency, 

cps 

Results using 
influence coefficients 
of table 1(b) 

Results using 
influence coefficients 
of table III 

Frequency, 

cps 

Deviation, 

percent 

Frequency, 

cps 

Deviation, 

percent 

First 

108 

107.9 

-0.09 

108.7 

0.65 

Second 

253 

214 

-15.3 

243 

-3.8 

Third 

342 

361 

5.5 

370 

8.1 

Fourth 

491 



463 

-5.7 


From each set of frequency and deflection calculations on wing 1, 
the node line associated with each calculated frequency was obtained 
graphically from cross plots of the computed deflection of the wing in 
that mode. The computed node lines which are shown in figure 3 appear to 
be in good agreement with the measured node lines. 
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In figure 5, plots of the calculated vertical deflections of each 
wing strip of wing 1 in the first three coupled modes computed from the 
influence coefficients of table l(b) are presented (dashed lines). . The 
deflections shown have been normalized by the extrapolated deflection of 
the tip of the quarter-chord line. The deflection curves show marked 
chordwise bending in the second and third mode for some of the strips. 

In order to simplify the theoretical flutter-speed calculations the chord- 
wise bending of the strips was removed by fitting a root-mean-square 
straight line to each of the strip deflection curves (solid lines in 
fig. 5). Since strip 5 has only two influence coefficient stations and 
strip G has one influence coefficient station and an extrapolated slope , 
the calculated and linearized deflection curves for these strips in 
figure 5 coincide. 

The vertical translation of the strip reference axis and rotation 
about the reference axis for each wing strip as obtained from the linear- 
ized deflection curves are presented in table Il(c) for each of the first 
three coupled modes. The values in table Il(c) of the angular deflection 
about the strip reference axes are normalized by the extrapolated values 
of the angular deflection of the streamwise tip chord about the tip of 
the quarter-chord line. 


Calculation of Reference Flutter Speeds 

Theoretical or reference flutter speeds were calculated for wing 1. 
These reference flutter speeds were computed by the use of streamwise two- 
dimensional incompressible aerodynamic coefficients in a coupled modal 
analysis. The frequencies used in the analysis were the first three meas- 
ured coupled natural frequencies and the mode shapes used were those of 
table Il(c) which were obtained as discussed in the previous section. 

The influence coefficient stations on each strip lie along a stream- 
wise strip chord as shown in figure 4. Each strip reference axis is the 
line normal to this strip chord at its quarter-chord point. In the ref- 
erence flutter-speed analysis , since the strip reference axes were effec- 
tively the same as if the wing quarter-chord line had been chosen for the 
reference axis, there was a simplification in the aerodynamic terms. 

The equations of motion and the procedure for computing the coeffi- 
cients of the flutter stability determinant are given in the appendix. 
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RESULTS 


General Comments 

For each of the flutter points only the modified wing panel fluttered. 
The other panel of solid 2024 aluminum was tested simultaneously, hut it 
was too stiff to flutter in the density range of the tunnel. 

The operating characteristics of the tunnel were such that frequently 
during a single test run (a test run is defined as one operation of the 
tunnel from valve opening to valve closing) the tunnel operating curve of 
dynamic pressure as a function of Mach number intersected the wing 
flutter-houndary curve of dynamic pressure required for flutter against 
Mach number more than once. In such instances, each point of intersection 
is presented in the data. 

Slightly more than half of the start-of-flutter points were readily 
determined from the oscillograph records. Each of these starts of flutter 
was characterized hy a change from random wing motion to continuous sinus- 
oidal oscillations accompanied hy an increase in oscillation amplitude. 

For these flutter points, when hoth sets of strain gages were operating, 
the frequencies of hoth strain gages were the same at the start of flutter. 
For the remainder of the flutter points a period of intermittent sinus- 
oidal type of oscillation preceded continuous flutter and obscured the 
exact start of flutter. Such periods are designated low-damping regions 
as in reference 1 inasmuch as the sum of the aerodynamic and structural 
damping is near zero. Where low damping occurred, two data points were 
selected: one point near the start of the low-damping region and the 

other near the start of continuous flutter following the low-damping 
region. Both data points are presented in the tables and figures. 


Presentation of Data 

The results of the investigation are listed in table IV. The first 
four columns of the table contain a description of the chronological 
behavior of the wing during each test run. The first column gives the 
wing identification number, the second column the number of the run, 
and the third column the chronological number of each data point during 
each run. The fourth column contains code letters (defined in table IV) 
which describe the behavior of the test wing panel at the time of each 
data point. 

Some of the experimental results tabulated in table IV are plotted 
as functions of Mach number in figures 6, 7, and 8 for both wings 1 
and 2. Figure 6 is a plot of dynamic pressure at flutter; figure 7 is 

— — — j and figure 8 is a plot of the flutter 


a plot of the parameter 
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frequency normalized by the third natural frequency, which is designated 
as the predominantly torsion frequency. It should be noted that for the 
point at Mach number 1.28 in figure 7 the corresponding ratio of frequen- 
cies has not been plotted in figure 8 because the flutter frequency was 
not obtained. A typical history of Mach number and dynamic pressure 
during a tunnel run is shown for wing 1 in figure 9 • Experimental results 
normalized by analytical results are shown for wing 1 in figures 10 
and 11. Figure 10 shows the variation of flutter-speed ratio with Mach 
number and figure 11 depicts the corresponding variation of the ratio of 
experimental to calculated flutter frequency with Mach number. In 
figures 6, 7, and 10 the low-damping regions are indicated by dashed 
lines which extend from the start-of -low-damping point (marked only by 
the lower end of the dashed line) to the continuous -flutter point (marked 
by a symbol at the upper end of the dashed line ) . 


DISCUSSION 


As shown in figure 6, flutter was obtained on wing 1 at Mach numbers 
from 0.79 to 1.28 and values of dynamic pressure from 878 to 4,296 lb/sq ft. 
Wing 2 was stiffer and slightly heavier than wing 1, so that the dynamic 
pressure required to flutter wing 2 was nearer the upper limit of the 
tunnel dynamic-pressure range. Consequently, only four flutter points at 
Mach numbers from 0.84 to 0.97 were obtained on wing 2 although it was 
attempted without success to flutter wing 2 at Mach numbers up to 1.2 
with dynamic pressure up to about 4,500 lb/sq ft. The dynamic pressure 
required to flutter wing 2 at a given Mach number was about double that 
required for wing 1. (See fig. 6.) Table IV shows that at a given Mach 
number the density and hence the flutter mass ratio differed for the 

two wings by a factor up to 2. (Compare, for instance, the values of 
wing 1, run 13, point 2 with those of wing 2, run 2, point 2 in 
table IV. ) It can also be seen from figure 6 that for wing 1 at low 
supersonic Mach numbers there is a very sharp rise in dynamic pressure 
required for flutter. These data can be discussed with more facility 
by making use of figure 7. In figure 7, for the two wings investigated, 
the data obtained are shown to be correlated by the parameter 

, which was also employed in reference 1. This fact indicates 

V e 

that, at a given Mach number, the flutter-speed coefficient 

for the two wings varies nearly linearly with the square root of the 
mass ratio for these tests. Figure 7 shows for wing 1 at a Mach number 

of about 1.05 a sharp rise in — Ve -- which corresponds to the similar 

b rVe 


Ve 
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rise in in figure 6. This sharp increase is interpreted as a change 

in flutter mode. Additional evidence that this shift represents a change 
in flutter mode is presented in figure 8, which shows a distinct jump in 
the ratio of flutter frequency to torsion frequency at the same low super- 
sonic Mach number as the shifts in figures 6 and r J. Similar shifts in 
flutter mode on a delta-wing model are reported in reference 2. 

The possibility exists that a second shift in flutter mode on wing 1 
may have occurred at a Mach number of 1.24 as evidenced by a large jump 
in dynamic pressure at flutter (see fig. 6) and frequency of flutter 
(see fig. 8). The frequency ratio for the point at Mach number 1.28 was 
not shown since the frequency was not obtained from the record; however, 
the dynamic pressure at flutter for this point is given (fig. 6). The 
data available are insufficient to draw a conclusion and the faired 
flutter boundaries of figures 7 and 10 are drawn as though a second shift 
in flutter mode did not occur. 

The values of the flutter frequencies of wing 1 in the low-frequency 
flutter mode (175 cps to 204 cps) are between the measured first and second 
natural coupled frequencies, whereas the frequencies of the high-frequency 
flutter mode (575 cps to 417 cps) are, with one exception, between the 
measured third and fourth coupled natural frequencies. The one exception - 
wing 1, run 10, point 4 - shows a flutter frequency of 500 cps. 

Special note should be taken of the regime below the flutter boundary 
for the high-frequency flutter mode. During the tests, because of the 
tunnel operating characteristics, a start and stop of flutter in the low- 
frequency mode was always obtained before flutter in the high-frequency 
mode was encountered. This sequence of wing behavior is illustrated in 
figure 9 which is a history of tunnel dynamic pressure and Mach number 
during tunnel run 11 on wing 1. The circled points in figure 9 are taken 
from the data given for run 11 in table IV. On every run, as dynamic 
pressure increased in the regime between the stop of flutter in the low- 
frequency mode and the start of flutter in the high-frequency mode, some 
random and intermittent oscillations were noted which were followed by 
a significant intermittent response of the wing to tunnel turbulence 
(fig. 9). These intermittent oscillations of the wing were of approxi- 
mately the same frequency as the low-frequency flutter mode. As dynamic 
pressure increased, intermittent oscillations of the wing at a frequency 
near that of the high-frequency flutter mode superimposed on the low- 
frequency oscillations and gradually replaced them (fig. 9)* After the 
high-frequency intermittent oscillations, flutter in the high-frequency 
mode began and continued up to the highest value of tunnel dynamic pres- 
sure reached during the run. 

The reference flutter speeds calculated for wing 1 as described in 
the "Analysis" section had several limitations: The use of two-dimensional 
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aerodynamic coefficients was unrealistic on such a. low-aspect -ratio, 
highly swept delta wing. The aerodynamic coefficients were also 
incompressible ones which precluded the possibility of predicting the 
change in flutter mode which was probably caused by changes in aerodynamic 
loading with Mach number. In addition, the effects of chordwise bending 
on the aerodynamic forces, as well as on the inertia and elastic forces, 
were neglected. In this connection, it might be pointed out that the 
neglect of chordwise bending distorted the true mode shapes and destroyed 
the orthogonality relationship so that the off-diagonal terms of the mass 
matrix which were considered to be zero were actually not zero. Only the 
first three coupled modes were used in the analysis. Since the frequency 
of wing flutter in the high-frequency mode was, with one exception, between 
the frequencies of the third and fourth coupled modes it might be expected 
that, in addition to the use of correct aerodynamic coefficients, the 
inclusion of chordwise bending and the fourth and perhaps higher coupled 
modes in the analysis would be required to predict the flutter character- 
istics. The inclusion of chordwise bending might also have improved the 
correlation at low Mach numbers. 

In figure 10 the flutter-speed ratios for wing 1, which were obtained 
by dividing the experimental flutter speeds by the reference flutter 
speeds, have an average of about 0.72 at Mach number 0.80. The flutter- 
speed ratios increase with Mach number to a value of approximately 0.90 
at Mach number 1.05 where an abrupt shift in V e /Vp, to a value of about 

1.15 occurs. After the shift, which is attributed to the change in 
flutter mode, the flutter-speed ratios increase steadily with Mach number 
to a value of about 1.35 at Mach number 1.28. As shown by the flutter- 
speed ratios of figure 10, the reference flutter speeds were too high 
(unconservative) for the low-frequency flutter mode by 11 to 28 percent 
and too low for the high-frequency flutter mode by 11 to 35 percent. 

Figure 11 shows that agreement between experimental and calculated 
flutter frequencies in the low-frequency flutter mode (values of cOg/cr^ 

around 0.95) is somewhat better than the agreement between experimental 
and calculated flutter speeds. However, in the high-frequency flutter 
mode the calculated flutter frequencies were much too low as shown by 
values of <x> e /co^ of about 1.9* 


CONCLUSIONS 

From a transonic flutter investigation of two 64° sweptback delta 
wings having streamwise ribs and orthogonal, spanwise spars but different 
stiffnesses and masses, the following conclusions were obtained: 
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1. At a given Mach number the values of the mass ratio for the two 
wings differed by a factor up to 2, but the data were correlated by use 
of the parameter consisting of the flutter-speed coefficient divided by 
the square root of the mass ratio. 

2. On the more flexible wing both the dynamic pressure required for 

flutter and the flutter frequency suddenly increased by a factor of about 
2 at a Mach number of approximately 1.05: This increase is interpreted 

as a change in flutter mode. 

5. Reference flutter speed calculations made for the more flexible 
wing using streamwise two-dimensional incompressible aerodynamic coef- 
ficients in a coupled modal analysis were too high for the low-frequency 
flutter mode by 11 to 28 percent and too low for the high-frequency mode 
by 11 to 55 percent. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , September 15, 1956. 
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APPENDIX 

EQUATIONS OF MOTION FOR REFERENCE FLUTTER-SPEED CALCULATIONS 


The equations of motion used for the reference flutter-speed calcu- 
lations are derived in reference 3 and are as follows: 


1 A i 


.2 . 


1 - (2) U - * 1 ) 


+ c llf &l + c 12^2 + C 13^3 = 0 



f 




C 21^1 + | 


[-(*)< 

4 + %) 

4“ C > 

-> 




r~ 

x 2 


c 3lll + C 52 l2 + ‘ 

r 

> 

V>J 


f) 4 + 

+ c 55 ► 


l 2 ^23^3 - ® 


b = ° 


where f-.(t) = | • e lGSb is the generalized coordinate which is a function 

1 1 y O 

of time, the amplitude of which expresses how much of each normal mode is 
included in the general vibratory motion. 


The necessary and sufficient condition that solutions for the simul- 
taneous equations of motion exist (other than = Ip = ^3 = °) is that 

the determinant of the coefficients equals zero. 

The coefficients in the equations of motion may be broken down as 
follows : 


Ai = Ti + 4*2 + T 3 
A 2 = 94 + T 5 + T 6 
A 3 = T 7 + 98 + T 9 
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where the 

values of cp for wing 1 

are as follows: 

*1 = y 

mh^b = 0 . 00011^99 

cp 5 = £ = 0.00000024200 

U'nI 

11 

OJ 

& 

I a o , 1 2 6 = 0.0000018522 

cp 6 = 2 £ S^gOgb = -0.0000027754 

95 = 2 

y s a h l a 'l & = 0.000012114 

cpy = ^ mhj 2 & = O.OOOO 59835 

% = Y 

mhg 2 6 = 0.00010864 

9 8 = Y I a a 5 2& = 0.000015351 



99 = 2 ^ S a h 3 a 3 6 = -0.000026389 


and the summations are over wing strips 1 through 6. The terms 

are: 
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C33 = JtP 
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~2 $ 31 + *32 
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where F r and G r are the Theodorsen functions corresponding to t=— 


by.0) 

and k-r = -1— and the values of 
r V 

*1 = ^10 + ^12 + 8 ^15 

= 0.000^9930 

<S> 2 = 2b r 2 cp llf = 0.00056098 

0 3 = 2b r (cp 11 + 9 13 ) = 0.0021970 

\ = h r ^cp 13 + 9p6) = 0.00015015 

*5 “ ’’17 + 2 *19 + 2 *22 + 5 *23 

= -0.000335^8 
<3> 6 = 2b r 2 cp 2 i = 0.00032259 
$7 = 2b r (q> 1 g + cp 2 o) = -0.00062021 


for wing 1 are as follows: 

$9 = ^25 + \ 9 27 + 2 9 30 + | 9 31 
= o.oooi8o4o 


® = 2b r cp 29 = -0.00055175 

*11 ' 2 M* 2 6 + *28 ' 

) = 0.00062875 

*12 - b r(*28 + *32) 

= -0.00025036 

<E> 13 = 2b r 2 9 5 ^ = -0.00025457 

*14 - 2b r(*L8 + *33' 

) = -0.00094354 

*15 - b r(*33 + *24) 

= -0.000094898 


$ 8 = b r (^20 + ^4 ) = °* oooo66 t 68 ® l6 = 935 + 937 + | 94 o = 0.00049528 
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$ 17 = 2b r 2 cp 59 = -0.000027012 

$23 — ^ > r(^ > 26 ^50) — 0.0011821 

*18 = 2b r(^6 + ^38 ) = 0.0014210 

$ 26 = b r( c P50 + ^32) = 0.000026320 

$ 19 = b r (^38 + ^41) = -0.000022526 

0> 27 = 2b r 2 cp 5 3 = 0.00017780 

*20 = %= + § V + 2 ^7 + 5 <p 48 

$ 28 = 2b r (cp^ + 952) = -0.00052604 


= -0.00012882 


*2i = 2b r 2 cp 46 = 0.00078900 

$ 29 = h r ^2 + ^49) = 0.000017796 

\ 2 = 2b r( c P43 + ^45) = 0.000175 66 

^30 = ^54 + ^56 + 8 ^59 = 0.00012599 

$23 = b r (^45 + ^49) = 0.00026864 

$ 51 = 2b r 2 cp 5 Q = -0.00036476 

$24 = 2b r 2 cp 51 = 0.00029982 

$ 52 = 2b r (cp 55 + 9 57 ) = 0.00044084 
$33 = b r (9 57 + 9 6o ) = -0.000020573 


The following values of cp for wing 1 can be evaluated from table II 
(the summations are over strips 1 through 6) : 


cp 1Q = £ b 2 h x 2 6 = 0.00042157 

9 i 4 = bh 1 a 1 5 = 0.0045975 

cpn = bh l 2& = 0.0039323 

q>l 5 = ^ b 4 a 1 2 6 = 0.000016207 

cp-^2 = Y b^h-jO^b = 0.000071649 

cp l6 = b 5 a 1 2 S = 0.000092856 

cp 15 = £ b 2 h 1 a 1 b = 0.00051504 

9 17 = b^l^b = -0.00031040 
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9is - 

bh-^hgb = 

-0.0014930 

*33 * I 

b 2 h 2 a x b = -0.00041701 

^19 = 

'y b^h-j^c^S 

= 0.000026839 

*34 - ] 

, bh 2 a- L b = - 0.0020863 

II 

O 

OJ 

a- 

X b 2 h 1 a 2 6 

= 0.00023751 

*55“ ] 

7 b 2 h 2 2 b = 0.00046584 

cp 2 l = 

'y bh^a^b = 

: 0.0026438 

936 = ^ 

> bh 2 2 b = 0.0029824 

cp 22 = 

'y b^h 2 a -|_6 

= -0.000080477 

*37= ] 

) b 5 h 2 a 2 b = - 0.000021208 

9 2 3 = 

b^a^agb 

= o.ooooo4634o 

*38= ; 

y b 2 b 2 a 2 b = -0.00010593 

cp 24 = 

X b^a^b 

= 0.000032808 

*39 = : 

^ bh 2 a 2 b = - 0.00022138 

^25 = 

^b 2 h x h 5 b 

= 0.00024453 

*40 = 

Vb^a 2 2 b = 0.0000017402 

*26 “ 

^ bh-jh^b = 

= 0.0020758 

*41 - 

^b 5 02 2 & = 0.000014730 

9 27 = 

X b^h^b 

= - 0.00015097 

%2 = 

^b 2 h 2 h 5 b = -0.00025049 

*28 = 

X b 2 h 1 a^b 

= - 0.00080303 

943 = 

^bh 2 h 5 b = -0.00078501 

^ 29 = 

^bh^b = 

= - 0.0045219 

944 = 

£b 5 h 2 a 5 b = 0.00023529 

II 

O 

rO 

a- 

X b 5 h 3 a l & 

= 0.000056820 

945 = 

^b 2 h 2 a 5 b = 0.0011406 

CPji = 

X b4a 1 a 3 5 

= - 0.000045483 

946 = 

^bh 2 a 5 b = 0.0064663 

*32 = 

V b^a a b 

z_ 13 

= - 0.00021057 

947 = 

^b 5 h 3 a 2 b = 0.000015340 
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’’Ml - ] 

r b^agOjb = -0.0000097254 

\ 

r b 2 h 5 % 

; 

b^o^a^b = -0.000052972 

”55 - ] 

— 2 

> bhj b ■ 

Vj = , 

b 2 h c^b = 0.00031713 

= ] 

[ b Vj 1 

II 

1 1 

Lf\ 

e- 

bh^a-jb = 0.00 24572 

*>7‘ ] 

r b2h 3 a 3 i 

cp 52 = 

V b^c^b = 0.00012502 

^8 ■ ] 

> bh^cx^b 

953 = 

bh^ctgb = 0.0014572 

^59 = , 

^ b^ 2 b 



’’So “ 

V b^a/b 


= 0.00024409 

= 0.0016278 

- = -0.00017737 

' = -0.00073342 
= -0.0029894 
= 0.00015805 


0.00065213 
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TABLE I.- INFLUENCE COEFFICIENTS ON WING 1 



ro 
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TABLE II.- PHYSICAL PARAMETERS 
(a) Wing properties 


Parameter 

Wing 1 

Wing 2 

A, deg 

63.6 

64 

Span, ft 

1.092 

1.085 

Panel span, ft 

0.421 

0.4i8 

b r , ft 

0.247 

0.244 

X cg , avg 

0.442 

— 

8h> ay g 

0.0283 

0.0216 



0.3158 

0.2991 



0.7398 

O .7196 

Exposed panel mass, slugs 

O.OO 587 

0.00681 

fl> cps 

108 

160 

f£, cps 

255 

385 

cps 

342 

555 


(b) Measured mass properties of wing 1 


Strip 

mb, 

slugs 

slug-ft 2 

^aP> 

slug-ft 

5, ft 

b, ft 

a, ft 

1 

0.00184 

0.000187 

0.0003435 

0.056 

0.438 

0.1866 

2 

.00188 

.0001219 

.000268 

.0833 

.360 

.1434 

3 

.00120 

.0000430 

.000128 

.0833 

.276 

.1067 

4 

.000674 

.00001297 

.0000483 

.0833 

.192 

.0716 

5 

.000241 

.00000173 

.0000084 

.0833 

.108 

.0358 

6 

.0000413 

.00000005 

.0000009 

.0311 

.0521 

.0208 


(c) Computed deflection properties of wing 1 


Strip 

h l 

h 2 

h 5 

a i 


a 3 

1 

0.00021 

- 0.0017 

0.0044 

0.00098 

0.00015 

-0.0064 

2 

.011 

-.027 

.039 

.0267 

-.0019 

-.1250 

3 

.063 

-.142 

.156 

.1070 

.0194 

-.4984 

4 

.194 

-.284 

.158 

.2597 

.0972 

-.8955 

5 

.479 

-.242 

.116 

• 5970 

.5022 

-. 66 29 

6 

.862 

.596 

.662 

.9311 

.8892 

.5918 


i 
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TABI* m.- SUTURAL INFLUENCE COEFFICIENTS d™> AT CENTER OF GRAVITY OF SECTIONS OF WING 1, tt x 10 5 



ro 

V>J 
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TABLE IV.- COMPILATION OF TEST AND ANALYTICAL RESULTS 


Wing panel behavior code: 

F - flutter 

E - end of flutter (dynamic pressure increasing) 

D - low damping 

Subscripts : 

1, 2, 3 - associated with first, second, or third occurrence 
of flutter during test run 



WING l 


i 



Figure 1.- Sketch of 64 delta wings showing construction details and 

dimensions . 
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Figure 2.- Variation of ratio of average thickness (over flat part of 
wing) to chord along wing span for wings 1 and 2. 




Figure 3*- Measured and calculated coupled natural vibration frequencies 
and node lines, shaker location, and strain-gage positions. 


ro 
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NOTE; All dimensions in inches 


Strip 

5treamwise distance 
trailing edge to Strip 
reference axis, m. 

\ 

7.89 

2 

6.48 

3 

4.97 

4 

3.46 

5 

1.95 

6 

,94 


Mast 

wing 


es of 
sections 


Section 

Moss 
(slugs x 10 s ) 

1 

63.60 

2 

37.67 

3 

34.41 

4 

41.77 

5 

6.43 

6 

37.64 

7 

48.79 

8 

41.92 

9 

50.56 

10 

8.45 

11 

33.51 

12 

37,42 

13 

41,24 

14 

7.48 

15 

26.37 

16 

34,69 

17 

6,58 

18 

18. 14 

19 

5.65 

20 

4.13 



Wing reference ans 
( quarter - chord line) 

-Strip reference axis 


• - Points at which influence 

coefficients were measured 

Points at which influence 

coefficients were interpolated 

□ - Measured eg position of 
Section 


Numoers locate influence 
coefficients and associated 
sections of wing 


Typicol section 
boundaries (sect. 16) 


~,206 


4.- Sketch of wing 1 showing strips, influence -coefficient sections, 
stations, and section center-of-gravity positions. 


Figure 


1 


__y 

y (tip of c/U line) 



(a) First mode - calculated frequency 107*9 cps. 


<8 


Figure 5* _ Calculated strip mode shapes and root -mean- square strip mode 

shapes of wing 1. 
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y (tip of c/k line) 
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(b) Second mode - calculated frequency 214.4 cps. 


Figure 5*- Continued. 
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(c) Third mode - calculated frequency 3^0. 9 cps. 


Figure 5 


Concluded 
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Figure 6.- Variation of dynamic pressure at flutter with Mach number for 

wings 1 and 2. 



Figure 7 .- Variation for wings 1 and 2 of the parameter 

^r^aFe 


with Mach number. 
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Figure 8.- Variation of ratio of flutter frequency to measured predomi- 
nantly torsion frequency with Mach number for wings 1 and 2. 
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Figure 9-~ Variation of tunnel dynamic pressure with Mach number for a 
typical tunnel run on wing 1, during which flutter apparently occurred 
in two modes. 
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Figure 10.- Variation of flutter speed ratio with Mach number for wing 1. 
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Figure 11.- Variation of ratio of experimental to calculated flutter 
frequency with Mach number for wing 1. 
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